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We study the large-Q^ behavior of the meson form factor Fm{Q^) constructed using the holo- 
graphic hght-front wave functions proposed recently by Brodsky and de Teramond. We show that 
this model can be also obtained within the Migdal's regularization approach ("meromorphization"), 
if one applies it to 3-point function for scalar currents made of scalar quarks. We found that the 
asymptotic behavior of Fm{Q^) is generated by soft Feynman mechanism rather than by large 

transverse momentum dynamics, which causes very late onset of the asymptotic behavior. It 

becomes visible only for unaccessible momenta > 10 GeV^. Using meromorphization for spin- 
1/2 quarks, we demonstrated that resulting form factor has l/Q* asymptotic behavior. 

Now, owing to the late onset of this asymptotic pattern, F^"^"" (Q^) imitates the behavior in 

the few GeV^ region. We discuss analogy between meromorphization and local quark-hadron duality 
model for the pion form factor, and show that adding the 0{as) correction to the spectral function 
brings in the hard pQCD contribution that has the dimensional counting behavior at large Q^. 

At accessible Q^, the 0(as) term is a rather small fraction of the total result. In this scenario, the 
"observed" quark counting rules for hadronic form factors is an approximate phenomenon resulting 
from Feynman mechanism in its preasymptotic regime. 
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I. INTRODUCTION 



Experimental evidence that form factors of hadrons consisting of n quarks have behavior close to 1/{Q'^)'^~^, 
provokes expectations that there is a fundamental and/or easily visible reason for such a phenomenon, scale 
invariance being the most natural suspect. Indeed, for an elementary fermion, EM form factor is constant, 
(P'|J'^|P) = u{P')^^u{P). When (n — 1) spectator quark fields with dimension (mass)'^/^ are added to the ini- 
tial state \P) and the final state the extra kinematical factors ~ [u(P)w(P')]"~^ take only worth of (mass)'^"~^\ 
and a constant of dimension (mass)^'^"^^^ is needed to take care of the rest. If, except for this overall constant, no other 
dimensionful constants can show up for large , then form factor has the quark counting rule l/((5^)"^^ behavior 
A specific dynamical mechanism ^^i] that produces a scale invariant behavior is provided by hard rescattering in 
a theory with dimensionless couplinig constant. After advent of QCD, pion and nucleon form factors were calculated 
within the hard scenario [1, 0, [1], and it was realized that perturbative QCD predicts, in fact, the (as/Q^)"~^ 
asymptotic behavior. In the pion case, the prediction is P^rlQ^) = (2as/7r)so/Q^, where sq = 47r^/^ ~ 0.7GeV^ is 
a constant close to « 0.6 GeV^. This indicates that the pQCD asymptotics is not the large-Q^ limit of the phe- 
nomenologically successfuU VMD fit P^^^(Q^) ~ 1/(1 -I- Q"^ /m^p), but rather looks like 0{as) correction to it. Also, 
the smallness of as/'"' undermines attempts to describe available data solely by pQCD hard mechanism. During the 
last years, the growing consensus is that at available , form factors are dominated by soft contributions described 
by nonforward parton densities J-{x, Q'^) Q (or generalized parton distributions iJ(x,^; Q^) for zero skewness 5), and 
successful fits were obtained [1, H, [l3| in models with T{x,Q'^) = f{x) e^'^'^^^^ having exponential behavior for large 
at fixed x. If g(x) vanishes for x 1 like (1 — a;)", a powerlike asymptotics F{Q^) ~ (1/0^)^''+^)/" in this case 
appears only after integration over x, i.e., it is governed by the so-called Feynman mechanism |11| . and is determined 
by the a; ^ 1 behavior f{x) ^ (1 — x)^ of the preexponential function f{x), in contrast to the hard mechanism for 
which the subprocess amplitude already has the (l/Q^)"""'^ power behavior not affected by subsequent x-integrations. 

Another new development is related to applications of AdS/CFT construction to QCD and claims [H, 111] that 
this framework (often referred as "AdS/QCD") provides a nonperturbative explanation of quark counting rules. In 
this scenario, they reflect the conformal invariance of the 5-dimensional theory, in particular, the power behavior of 
the normalizable modes $(C) at small values of the 5th coordinate C. Parametrically, the prediction has the form 
(A^/Q^)"^^, without any accompanying a"^^ factors. Given explicit expressions for ^(C), with the value of A fixed 
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from fitting the hadron masses, it is straightforward to check the structure of AdS/QCD results for form factors and 
their potential to describe the features of existing data. This is one of the goals of the present investigation. Another 
is to study the recently proposed interpretation of AdS/QCD results in terms of light-front wave functions, 
which opens a possibility to find out whether, in terms of the light-cone momenta x,h±, the AdS/QCD quark 
counting corresponds to large-k_L hard mechanism or, as we will show, to the x — > 1 soft Feynman/Drell-Yan pjl [isj 
mechanism. In view of yet another recent observation [3| that some of the results of the AdS/QCD approach coincide 
with those of Migdal's program llTl (that starts with a perturbative correlator and substitutes its cuts by hadron 
poles), we apply the extension |18j of this "meromorphization" idea to the 3-current correlators, and establish a 
connection between this approach and holographic light-front wave functions of Ref. [l^ . 

Finally, we discuss analogy between the meromorphization procedure and the "local quark-hadron duality" model 
[Tgl that succesfuUy describes the pion form factor data and gives a unified description of its soft and hard parts. 
The soft term F^^{Q'^) in this approach dominates at accessible Q^, but has the 1/Q^ asymptotic behavior. Due to 
its late onset, the curve Q^F^^{Q'^) has a wide plateau in a few GeV^ region, i.e., F^^{Q'^) imitates there the 
behavior. Thus, the desired 1/Q^ result for accessible is obtained because the nonperturbative term has a faster, 
l/Q^ asymptotic fall-off. On the other hand, the curve Fm (Q^) for the meson form factor in the holographic model 
of Ref. 14] monotonically increases to its asymptotic value ~ 2.6 Wp, and is far from being flat in the whole accessible 
region < lOGeV^. 



II. HOLOGRAPHIC WAVE FUNCTIONS 



We will need some elements of the derivation of the holographic wave functions proposed in Ref. In the 

hard- wall approximation (T^ , the expression for the elastic form factor in the holographic duality model is given by 

P(Q') = J^ ^*P'(C)J(Q,C)<i>p(C) , (1) 

where J{Q, C) — CQKi{CQ) is the "nonnormalizable" mode describing the EM current, and <I>p(C), <I>p'(C) come from 
the "normalizable" modes describing initial and final states, with = CC,"^ JLi^L.kC^) for a qq state with orbital 
momentum L = 0, 1, . . . and radial number = 1, 2, . . ., and /3l,/c being the kt\\ root of the Bessel function Jl{x). On 
the other hand, in the light-cone (LC) formalism, 

F(Q2)^ f dx f d^rj^e'^^-^^B{x,Tj^)= f dxJ^{x,Q^), (2) 
Jo J Jo 

where B{x,ri±) is a parton density function 20] that depends on x, the light-cone momentum fraction of the active 
quark, and ri±, the Xj-weighted transverse position of spectators. The function B{x,r]±) accumulates information 
from all Fock components and its ri± Fourier transform gives the generalized parton distribution (GPD) ^{x, Q^). In 
particular, the 2-body part of a meson form factor is given by 



Fi2){Q^)^ dx J d^b^e^b^q^ ^2{x,h^) , (3) 
where x = 1 — x. With the wave function depending on b_L through b = \hj_\ one obtains 

F(2){Q'') = 2^ dx bdbJoixbQ)\^2ix,b)\'' = 2n ^ zdzJo(^'^zQ^\^{x,z)W (4) 
where z = \lxxb and wave function was written as ^2(2:, ^) — <t>{x, z). Noticing that the integral 



^ dxJol^-^^-^zQj = zQKi{zQ) = )CiizQ)^ J{Q,z) (5) 
gives J{Q,z), and assuming that )(/)(a;, z)]^ = xxx^(z), LC formula (|4]) is cast into the form 

/•OC 

F(Q2) ^271 zdz /Ci(zQ) x\z) (6) 
Jo 
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that converts into Eq. ([T]) if one identifies $(C) C^x(C)- I'^ general case, z is introduced through \r]^\ = z^/x/x, 
and one assumes that B{x,'r]±) = {x/x)x^{z). This gives the "holographic" model ^14] for light-front wave functions. 
By construction, they should guarantee the Fm{Q^ = 0) = 1 constraint, i.e., they are effective wave functions (see, 
e-g-, mi) normalized by 



dx I d^h±_\^cs{x,hA 



(7) 



unlike the two-body wave functions ^'2(2;, bj^) normalized to the probability of finding the meson in the qq Fock state. 
We will focus on the lowest meson having L = Q^k — \ and mass M = /3o.iA. Then 



where /3o,i Ji(/3o.i) ~ 1-2. The kx counterpart of this wave function is given by 

T / , N 1 /" ,9, ih W T / , X M JniBn ik±/\/xEM) 



(8) 



(9) 




FIG. 1: Momentum wave function My/nxx'i' m{x, kx) (lower curve, blue online) as a function of k\/M'^xx; for comparison 
also shown e'K\>{—k\/2M'^xx) (upper curve, red online). 

Since /3o,i is a root of Jq{z), the momentum wave function ^^(a;, k^) has no singularities for k'^/xx = M^, and it 
has zeros when k\/xx coincides with squared masses of higher states. For large k±, wave function oscillates with the 

magnitude decreasing as 



^'A/(a;,kj 



M y^-^ cos(fc^/\/xxA - 7r/4) 



TTkl 



\J k^/\/xxA 



(10) 



This result contradicts the statement made in Ref. that the AdS/QCD construction corresponds to a purely 
power-law large-fc_L behavior 'ip{k±) ^ (1/fc^)"^^ for the wave function of n-particlc bound state. 



III. FEYNMAN MECHANISM 



We can get the large-Q^ asymptotics of the lowest state form factor 



Fm{Q^) 



2A'P 



Q/A 



Q^Po.iJiiPo.iW Jo 

by Taylor expanding Jq{^M/Q) and neglecting 0{e^'^^^) terms from the upper limit of ^ integration: 



Q2[/3o,iJi(/3o,i)]' 



4Af2 9 /4M2 



0{e 



0.64 



1-HQV4M2 



(11) 



(12) 
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Though this result has a monopole-hke structure, the scale 4M^ is evidently too large. With M = nip, the curve 
for Q^Fm{Q^) is far from being flat in the accessible region < 10 GeV^. Note also, that though the large-Q^ 
behavior of Fm{Q^) is determined by the small-^ behavior of <&(C), when C is interpreted as \fx%h^ the value C = 
may correspond to a; = 1 rather than to large transverse momenta. Let us see which mechanism is responsible for the 
l/Q^ asymptotics in terms of the x, kj^ Drell-Yan formula [15i] 



dx / d'^k±'^li{x,'k±+xq±)'^Mix,'kj_) 



(13) 



Q = |qj_|. With 5'Af(x,kj^) decreasing at large k_L, there are two possibilities [23] (see also |23|): 
a) finite x and small |kx|, e.g., the region |kj^| ^ a;|q±|, where m{x,^±) is maximal. This gives 



Fm{Q^) 



dx / rf^k^ |«'^,^(a;,a;q_L)*M(a;,k^)| + {* 4^ **} - 2 / dx {'^liix, xq±) ip{x 



(14) 



In the hard scenario, when the x integral is not dominated by the a; ^ 1 region, the large-Q^ behavior of the form 
factor repeats the large-fc_L behavior of the momentum wave function. The latter does not behave as 1/A:^: it oscillates 
with magnitude decreasing as l/k'^^, see Eq. (fTO]) . so we need to turn to the second possibility: 

b) X is close to 1, so that |xqj_| ~ |k^|, and |kj^| is small. Then both ^m(2^, k_L) and ^'^^(a^j k_L+xq^) are maximal. In 
Ref. (isj . it was argued that the dominant contribution comes from a;|qj^| ^ ni — const, and the large-Q^ behavior of 
the form factor in this scenario reflects the phase space available for such configurations. To make a specific estimate, 
we represent the form factor as the a:-integral ^ of GPD 



Po.iJi iPo,i) Jo 



(15) 



which in this case is a function G{o') of cr = xQ'^/xA'P (see Fig. [5]). The form factor is then given by 



FMiQ^) 



da 



(16) 



and the leading term is proportional to zeroth (T-moment of ^?(c), which does not vanish though the function 

G{cr) oscillates for large values of a. In fact, it is very small in the region above its first zero (at a ~ 9) while in the 
region below this zero Gic) is very close to the exponential e""^/^'^, the integration of which gives the correct coefficient 
for the 0{AP /Q^) term. Thus, the leading large-Q^ term is given by integration over cr < 10 or, returning to the 
x-variable, by the region x < lOM^/Q^ (we see again that asymptotic estimates can be used only for ^ 4M^). The 
X < 0{AP /Q^) result parametrically differs from the estimate given in Ref. [15]. The reason is that the transverse 
momentum k_L enters into the wave function through the combination ]kj^]^/a;x, and the a;-size of the dominant region 
is determined by ]kj^ + xq±\'^ /xxM'^ < const or \xq±\'^/xx < const Af^, which gives x < const M^/Q^. 




2 4 6 8 10 12 14 



FIG. 2: Generalized parton distribution Tm{x, Q^) (lower curve, blue online) as a function of xQ^/xM^ = a. For comparison, 
we also show exp(— o-/3) (upper curve, red online). 

Thus, we found that the large-Q^ asymptotics of the meson form factor in the model of Ref. [ij] is governed by the 
soft Feynman mechanism, with the power-law asymptotics determined by the x — *■ 1 behavior of the prefactor f(x) 
accompanying a decreasing function of xQ^/xA^. In the present case, f(x) = 1, which gives Fm{Q'^) ^ l/Q^- With 
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extra factor, the outcome of the x-integration would be (A^/(5^)^+^. 
but the parton distribution function: 



Clearly, the prefactor f{x) is nothing else 



/(x) = ^(a;, Q2 ^ 0) 



(17) 



i.e., the model of Ref. [l4| gives a constant, x-independent parton distribution f{x) — 1. 

If we approximate ^Mix, k^) by a Gaussian function '^'cix, k^) ~ exp[—k'^/2M'^xx]/\/x^ (see Fig. [T]), the relevant 
GPD Tg{x, Q^) can be calculated analytically. It is instructive to perform such a calculation in the impact parameter 
space. Then ^'g(x, b) ^ y/x^ exp[—b'^M'^xx/2] (in fact, a Gaussian wave function for a ground state is obtained in the 
hologaphic model of Ref. [13] that gives a linear law for (mass)^ of excited states), and integrating over all (including 



small) values of bj_ in the form factor formula ^ gives J-g{x,Q'^) = e 



a function that, for any fixed x, 



vanishes exponentially at large Q . The power-law asymptotics is obtained only after integrating over x, which gives 



n=l 



4M2 



1 - 2 



AM' 



UP 



(18) 



a result that has the 0{M'^ /Q"^) large-Q^ behavior and structure similar to that of Eq. (fT2|) . The crucial role 
of integration over the a; ~ 1 region in getting the power-law behaviour is evident: if the integration is restricted 
to a; < Xo, the outcome vanishes exponentially (like exp[— a;oQ^/4a;o-^^]) for large Q^. Thus, the GPD !Fg{x,Q^) 
corresponding to the Gaussian wave function has the same J-{x, Q^) — f(x) e~'^ ^^^'^ structure as those considered in 
Refs. Mi,[ii. 



IV. MEROMORPHIZATION 

One may question the conjectures of the model of Ref. [l3|, e.g., the interpretation of the holographic variable C as 
a particular product \fx%h of light-cone variables. We are going to show that the picture similar to that of Ref. [IJI 
emerges also within the approach related to Migdal's prograi n [TtI of Fade approximating the correlators n(p^) of 
hadronic currents calculated in perturbation theory. Recently [To] , it was demonstrated that some of Migdal's results 
coincide with those of the holographic approach. Migdal's program involves "meromorphization" 



n(p') = - / ds ^ Um ip') = n(p^) - —— / ^^^^^ ds (19) 

that substitutes the original correlator n(p^) by a function IiM{p^) in which the cut of the original correlator for 
real positive is eliminated by the second term in Eq. (|19p . with zeros of Q{p^) at timelike p^ generating the poles 
interpreted as hadronic bound states. Explicit Fade construction in case of correlators of currents like the scalar 
current ipip of scalar fields, vector current V'7^V' of spin-1/2 quarks, etc., gives Q{p^) ^ Jq{Po,i\/^/^)j with M 
being the mass of the lowest state. In the deep spacelike region {p^ = —P^ ~oo), the difference between the 
original expression and the approximant then vanishes exponentially like e"^''"'^^/^^. The coupling constant of 
the lowest state is given by 



fl, = lin^(M'-p')UMp') = r^^^^ds . (20) 



In the lowest order, p[s) — pq9{s), with po = I/IGtt for j = (pip, and po = Nc/127: for vector and axial currents 
M7^(75)d of spin-1/2 quarks. This gives 



7r/3o,i Ji(/3o,i) X M^-s APo,iJi{PosW ^^^^ 

For form factors, it was suggested [l^l to "meromorphize" the 3-point hmction T{p^,p2,Q^). The lowest-order 
triangle diagram has only the double spectral density p(si, S2, Q")- Building the function 

T{Pi,P2,Q )^T{pj_,p2,Q )+ 2\^^ 2\ / '^^i / ds2 2T7 T\ ' (22) 



ir^QipDQipl) Jo 'Jo isi-plKs2-pl) 
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one removes the cuts in the pf and P2 channels substituting them by poles at the same locations as in n^vidj^). From 
this expression, one can extract the elastic form factor of the lowest state by using 



fiiFMiQ') = Hm ]im (pi - M'){pi~ M')T (pi pi Q') 



1 



7r2[Q'(Af2)]2 



dsi 



p(si,S2,Q')g(gl)Q(s2) 
(Sl - M2)(S2 - M2) 



(23) 



The spectral densities p(si,S2,<3^) can be calculated [2l|,[2^ using the Cutkosky rules and light-cone variables in 
the frame where the initial momentum pi has no transverse components pi — {p^ — V,Pi = si/V,0±}, while the 
momentum transfer q = p2 — Pi has no "plus" component, P2 = {V, (s2 + q^)/?-", q±}: 



P(S1,S2,Q^) 



'■1 , n{x) 
Pq I ax — — 



(fc_L + xq^)' 



XX 



(24) 



where x is the fraction of V carried by the active quark, and kj^ is its transverse momentum; pQ is the same as in Eq. 
()2ip . For the simplest correlator of three scalar ipip currents, the numerator factor is n{x) — 1/x. Taking i(p for 
the EM vertex gives n{x) — 1 since i xV^ . Then 



Jm^m W 



Pa 



XX 



Q(kl/xx) QiCk^+xq^)yxx) 



M2 - k5_/a;x ~ {k± + xq±)'^/xx 



(25) 



i.e., a LC form factor expression. Using Q{s) — Jo(/3o,i\/^/M) and Eq. (PT]) for , we obtain that the relevant 
wave function ^5vf'''''(a;, kj^) coincides with the momentum version ^ of the holographic wave function of Ref . . 
supporting interpretation in terms of light-cone variables x, b proposed there. 



SPINOR CASE 



Using spin-1/2 quarks and vector currents jajjp for hadronized vertices, one should deal with the amplitude 
T^p[pi,P2) and choose which tensor structure to consider. As a simple example, let us take the projection T^pn^vf^n^ , 
where n is a lightlike vector having only the minus component in the frame specified above. For p-meson form factors, 
this projection picks out the combination C(Q^) = Fi{Q'^)+kF2{Q^) — h^F3{Q'^), where k = Q"^ /2m?p. For comparison, 
pQCD expects that the leading l/Q^ behavior is provided by Fll{Q'^) = -Fi(Q^) - nF2iQ'^) + {k^ + 2k)F3{Q'^), with 
FiiQ^) ^ F2{Q^) ^ and F^iQ^) ^ (see, e.g., HI). Thus, -C{Q^) differs from Fll{Q^) only by Fi{Q^) 

and kF3{Q^) terms which are considered in pQCD as nonleading. The simplicity of this projection is that, for a 
spinor triangle, the numerator trace is given by the product of quark light-cone "plus" momenta xV, xV and xV . As 
a result, n{x) = 6xx. As we argued above, an extra x factor in GPD J^{x, Q^) should result in the l/Q^ behavior of 
the form factor at large Q^. To check this, we switch back to the impact parameter space representation and observe 
that the integral given by Eq. ^ changes into 



6 / dxxx Jo 
Jo 




(26) 



Hence, the form factor is given by an expression similar to Eq. ([TT|), but with /Ci(^) substituted by /C2(^). Both 
functions are equal to 1 for ^ = and behave as at large ^. Thus, one may expect that -F|/"^°'^(Q^) also behaves 
like l/Q^ for large Q^. But exphcit calculation gives 



Fir°'iQ' 



2M2 



Q2[/3o,iJi(/3o,i)]2 
2M2 



AP 



1 M2 ^ _ A .a O(MVQ^) 

2^ Q2 + Q4 576 ^ ne + ^^^'^ > 



0{e 



O-t-24 



288 



Af4 



2400 ■ 



o(m7q8) 



0{e 



(27) 



g4 Q6 

0, and the leading term has l/Q^ behavior. However, the next, 
7 GeV^, and the asymptotics establishes somewhere above 



The first term vanishes because d£^ JC2 (C) = 
correction exceeds it for all up to 12M2 
20GeV^. Now, this fact is very welcome: due to it, the i^^f '"""^(Q^) curve in the region of a few GeV^ imitates the 
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"power counting" l/Q^ behavior much more successfuhy than Eq. (fTTj) that displays its nominal asymptotics 
only well outside the few GeV^ region (see Fig. [3]). 

The miraculous cancellation of the moments for K2{£,) and K3{^) in ([27|l producing the Fl^"^°'{Q^) ^ M'^/Q'^ 
result can be traced to the 1/Q^ asymptotic behavior of the underlying double spectral density p{si, S2, Q^)- To this 
end, consider the double Borel transform of the three-point function [3) 111] 



^n,T2,Q^) = — dsi ds2 p{si,S2,Q^)e-''^'-''^' , (28) 

Jo Jo 

in which the power weights are substituted by the exponential ones. For the triangle diagram (Toj , 

^(ti,t-2,Q^) = ^ 3. r/ dxxxexp 

27r^(ri +T2) Jo 

It contains the xx factor (absent in the case of scalar quarks) that results in the behavior of $(ti,T2,Q^) and, 

hence, of the spectral density p{si,S2,Q^) = Nc 0{si) 9{s2) (si + S2)/2(5^ + . . . . 



XT1T2 

x{ti + T2) 



(29) 



VI. LOCAL QUARK-HADRON DUALITY 

Referring to the "established" behavior of meson form factors at large Q^, one primarily has in mind the data 

on the pion EM form factor which resemble the VMD monopole form Fj^iQ^) ^ 1/(1 + Q^/m^). In fact, the data 
are well described by our local quark-hadron duality model [19] that incorporates, in a simplified form, some ideas 
of Migdal's program and the QCD sum rule approach [27i] . From the latter, we borrow the observation that only 
the lowest state is narrow, and use the model spectrum: first resonance plus perturbative "continuum" starting from 
some effective scale s — sq. In other words, we transform the two-current correlator n(p^) into 



n(p^) ^ - / ds^-^^ n--(p^) = -— ^ + - / ^-^ds. (30) 

Then we try to reach the best possible agreement between the original and model correlators in the deep spacelike 
region of = —P^. For axial currents, Fm ~> /tt, M 'm-n ~ 0, pP'''''(s) — l/47r, and 

n(p2) - n^^ip^) ^Il + ^-Tds (31) 



p2 47r2 Jq s ~ p^ 

To eliminate the leading 1/p^ term in this difference, we should take = sqIAiv^ . As a result, the coupling and 
the value of sq are connected by the local duality relation 

r p^{s)ds= r pp''''\s)ds (32) 

Jo Jo 
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between the spectral density p7r(s) = Trf^5{s) of the lowest state and the perturbative spectral density pP*^''*(s), applied 
within the "duality interval" sq. Transforming the n""n^n^ projection of the 3-point function 



9 2 

PI Pi 



/ dsi 



1 - 0(si < So)(?(s2 < So) 



„pert 



(S1,S2,Q^) 



{si - pI){s2 - pI) ' 



(33) 



and requiring that the difference between TP''''*(p^,p2) Q^) and the model function vanishes faster than l/pfp^ for 
large spacelike Pi,P2, gives the local duality relation for the pion form factor [T^ 



Jo Jo 

Using Eq. ([M)) for /9P"''(si, S2, Q^) with n{x) — 6xx and = sq/Att'^ produces the DrcU-Yan formula 



Fi^°(g2) = — / dx f d^k^ e (ki < xxso) 9 ((kj 



a;q_Lj 



< xxso) 



(34) 



(35) 



with the effective "local duality wave function" \['^^(a:, kj^) = -^/B/ttso ^(^1 < xa;so) for the pion. In the impact 
parameter representation, "ii^^ {x,h±) = ^ydxx/ir Ji{by/xxso)/b. Taking the 6^0 limit gives the model distribution 
amplitude (p^^{x) = G/t^xx, that coincides with the asymptotic pion DA. Explicit expression for the pion form factor 
in the local duality model is known from [l9| : 



1 + Gsq/Q'^ 
(l + 4,so/Q2)3/2 



6sl 



210s4 



1008s^ 



(36) 



As expected, F^^((5^) behaves like l/Q"* for large Q^, but the expansion parameter of this series is again large 
~ 6so/Q^ ~ 4GeV^, so that the asymptotic l/Q'^ estimate should not be used at accessible Q^. However, in the 
region of a few Q^the F^^{Q'^) curve successfully imitates the 1/Q^ behavior and goes very close to existing 
and preliminary [29[ experimental data. 
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FIG. 4: "Holographized" local duality wave function ^7rso/6 *^°'''°'°(2;, kx) vs. k^/^/xE in GeV. 

The sharp cut-off form of 5'^'^(a;,k^) is a consequence of a simple model for the higher states' contribution. The 
resulting 6-space function \l/^°(a;, bj^) has a Bessel-type form, but goes beyond the first zero. "Holographizing" it by 
imposing the cut-off d{by/xxso < we found that its kj_ version has a smooth dependence on k'^/xx, qualitatively 

similar to that of Eq.([9]), with the lowest zero located at k\/xx « (1.26 GeV)^, i.e., "unexpectedly close" to the Ai 
position (see Fig. 2]). 



VII. HIGHER ORDERS AND TRANSITION TO PQCD 



Both in the meromorphization and local duality approaches one can include higher order as corrections to spectral 
densities. In particular, the two-loop calculation of p(si, S2, Q^) for axial and vector currents was performed in [33|. 
Among 0{ag) contributions, there are gluon-exchange diagrams whose asymptotic large-Q^ behavior is determined 
by the hard pQCD mechanism (see Fig. E]). As a result, the leading term of the spectral density corresponding 

to the T^pn^n°'n^ projection can be written in pQCD-like form, 

P.. {SUS2, Q') = 2.a. ^ dx dy 'jj^i^' + Oi^Q^) , (37) 
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where p{x, si) is the lowest order term for the x-unintegrated spectral density, 

N f / \ N 

P(^, 5i) = ^ y -5 (^51 - ^ j d^kix = ^ e{s,)xx (38) 

(and similarly for p(y, S2)). 




FIG. 5: Diagrams producing spectral density p{si, S2, Q^) at the one-loop and two-loop {0{as)) level. Outlined diagrams are 
responsible for the leading pQCD behavior. 

Integrating p{x, s) over the duality interval < s < sq and dividing result by nfj^ gives the local duality model 
(pY^{x) = dfyrxx = ip'^^ix) for the pion DA. Thus, substituting the 0{1/Q^) part of Pa^ (si, S2, Q^) into the local 
duality relation ((M)) gives the pQCD hard gluon exchange contribution i^P^*-^°((5^) = Sirasf^/Q^ calculated for the 
asymptotic shape of the pion DA. It is instructive to rewrite this result in the form FP^'-'^(Q^) = 2(so/Q^)(as/7r) 
that clearly reveals its nature as the as correction to the soft contribution ([36|l . with the {as/n) factor being the 
standard penalty for an extra loop. Using full expressions for Pa, (si, S2, Q^) one can get the local duality model 
prediction for F^^'"°'"^ at all Q^, see Ref. [s^. In fact, a very good approximation is given by a simple interpolation 
formula F^^'^°"-\q^) = (a^/7r)/(l -h QV2so) [IH between the = value F^^'^°"'\o) = as/n (that can be obtained, 
e.g., from the two-point result p{s) = S{s){l + as/ir) using the Ward identity) and the large-Q^ asymptotic behavior. 
With Os/tt « 0.1, the 0{as) term is a < 30% correction to the 0{a°) term in the < 4 GeV^ region (see Fig. [S]), 
and their sum is in good agreement with existing data. 
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2 4 6 8 10 

FIG. 6: Pion form factor in local quark-hadron duality model: Q^F^^'"'' (lower curve, red online), Q'^F^° (middle curve, 
green online), total contribution (upper curve, blue online). 



VIII. SUMMARY AND CONCLUSIONS 

In this letter, we studied the large-Q^ behavior of the meson form factor FmIQ"^) constructed using the holographic 
model of Ref. '14], and observed that, despite its asymptotic behavior, the combination Q^F]\i{Q^) is not flat 

in the accessible region < 10 GeV^. We also found that the asymptotic result is governed by the Feynman 

mechanism rather than by large transverse momentum dynamics. We discussed the meromorphization approach, in 
which the form factors are given by integrating the perturbative spectral density p{si, S2,Q^) with weights proportional 
to Jo{Po,i^/si/M)/ [si — M^), and showed that, if one uses scalar currents made of scalar quarks, the result coincides 
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with the expression given in Ref. [l4|. For spin-1/2 quarks, we studied a particular tensor projection of the 3-point 
function, and demonstrated that, due to an extra (1 — x) factor, the resulting form factor -F^f (Q^) has l/Q'^ 
asymptotic behavior. However, owing to the late onset of this asymptotic pattern, the combination 
rather flat in the phenomenologically important few GeV^ region, i.e., F^™°'' {Cf') imitates the 1/Q^ behavior. Then 
we presented the results for the pion form factor obtained in the local quark-hadron duality model, which corresponds 
to integrating the perturbative spectral density p(si,S2,Q^) with 9{si < sq) weights. Again, the lowest-order term 
has nominally the 1/Q* asymptotics, but it imitates the behavior in the few GeV^ region. Including the 0{as) 

correction term of p{si, S2, Q^) brings in the hard pQCD contribution having the dimensional counting behavior 
at large Q^. However, at accessible Q^, the 0(0;^) term is a rather small fraction of the total result, because of small 
as/n ~ 0.1 factor associated with each higher order correction. 

We did not discuss the nucleons form factors in the present paper, but we want to mention that the lowest-order 
perturbative density for spin-1/2 quarks is known [l^, and since the double Borel representation (see Eg. ipS)) ) in that 
case has (1 — x)^ factor, the resulting asymptotic behavior of p{si, S2,Q'^) is 1/Q^. But the local duality result for 
G\j{Q'^) closely follows the dipole shape of the data up to ~ 15 GeV^. 

Thus, we observe that the power of (1 — x) in these perturbative versions of the relevant parton densities f{x) 
increases with the number of quarks like fn{x) ^ (1 — a;)"~^, i.e., the probability that the total momentum of n — 1 
spectators is Xsp goes like a;"p~^, which looks quite natural. Due to Feynman mechanism, this formally gives (l/Q^)" 
behavior for form factors. However, because of the late onset of the asymptotic regime, the form factors imitate 
(1/Q^)"~^ behavior in a rather wide preasymptotic region. In this scenario, the "observed" quark counting rules is an 
approximate and transitional phenomenon dominated by nonpcrturbative, long-distance aspects of hadron dynamics. 
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